Quasihole dynamics as a detection tool for quantum Hall phases 
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Existing tecliniques for syntliesizing gauge fields are able to bring a two-dimensional cloud of 
harmonically trapped bosonic atoms into a regime where the occupied single-particle states are 
restricted to the lowest Landau level (LLL). Repulsive short-range interactions drive various tran- 
sitions from fully condensed into strongly correlated states. In these different phases we study 
the response of the system to quasihole excitations induced by a laser beam. We find that in the 
Laughlin state the quasihole performs a coherent constant rotation around the center, ensuring con- 
servation of angular momentum. This is distinct to any other regime with higher density, where the 
quasihole is found to decay. At a characteristic time, the decay process is reversed, and revivals of 
the quasihole can be observed in the density. Measuring the period and position of the revival can 
be used as a spectroscopic tool to identify the strongly correlated phases in systems with a finite 
number of atoms. 
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I. INTRODUCTION 

Strong correlations and anyonic excitations are the in- 
triguing properties of quantum states in two-dimensional 
systems exposed to strong magnetic fields. They show up 
in the context of fractional quantum Hall effect of elec- 
trons [11,0]. In recent years, the advances in techniques 
for cooling and controlling atoms have raised the hope 
that these interesting states might also be artificially gen- 
erated in systems of ultracold atoms 0]. This would al- 
low to experimentally confirm the fundamental theoreti- 
cal concept of fractional quantum statistics [J] , and open 
the door for topological quantum computation [5|. 

The key requirement for realizing such states is a 
strong external gauge field, which due to the electroneu- 
trality of the atoms has to be synthesized. Artificial 
gauge fields which are strong enough to bring the system 
into a regime, where only the lowest Landau level (LLL) 
is occupied, have already been generated by rotating a 
gas of *''Rb 0. The occurrence of strongly correlated 
states in the LLL regime then crucially depends on the 
ratio between trapping energy, favoring condensation in 
states with small angular momentum, and the strength 
of repulsive interactions, which tends to spread the atoms 
over a wide range of angular-momentum states. In a sys- 
tem of bosons interacting via a two-body contact poten- 
tial, this competition is known to restrict the Laughlin 
state [l| to a narrow region of parameters of extremely 
weak effective trapping, and thus close to the instabil- 
ity at the centrifugal limit ■ This drawback has so 
far hindered the experimental realization of the Laughlin 
state. It has led to the proposal of using laser-induced ge- 
ometric phases to mimic magnetic fields (cf. [13 )■ Such 
a method, experimentally proven in Ref. [T]| . allows for 
a precise tuning of the gauge field strength as required 
for reaching the Laughlin state. An experimental route 
to produce the Laughlin state could start with preparing 
the system in a condensate at zero angular momentum. 



L = 0. Then, stepwise transitions into states with higher 
angular momentum can be induced by adiabatically in- 
creasing the gauge field strength [l2| , until reaching the 
bosonic Laughlin state, characterized by L = N{N — 1) 
(in units of h) with N the particle number. 

An important question is then how to detect this 
state. Its zero compressibility or its constant bulk den- 
sity are characterizing features, but do not uniquely dis- 
tinguish the Laughlin state from other quantum liquid 
states. Moreover, in systems of only few particles these 
attributes may become quite unsharp, while experimen- 
tal progress in realizing Laughlin states of few particles 
has been reported (Tsj , and even small systems have been 
predicted to support bulk properties like fractional exci- 
tations Thus, looking for distinctive features, exper- 
imentally accessible even in small clouds, seems to be 
expedient. 

In this paper, we discuss a scheme for testing many- 
body quantum states in the LLL by piercing a quasi- 
hole into them. Experimentally, this can be achieved 
by focusing a laser beam onto the atomic cloud. After 
switching off this laser, the subsequent dynamics of the 
quasihole can be observed in the density of the system. 
We show that it yields relevant information about the 
underlying state. The defining property of the Laughlin 
state, being the densest state with zero interaction energy 
in a two-body contact potential, is found to be reflected 
in a decoherence-free dynamics of the quasihole. This 
is in clear contrast to the time evolution of a quasihole 
pierced into a state with L < N{N — 1). In this case, an 
interaction-induced dephasing delocalizes the excitation, 
visible in the density as a decay of the quasihole. We ex- 
plicitly consider a quasihole in the L = condensate, and 
in a Laughlin-type quasiparticle state. For these states 
we show, that the decay process is reversed at a charac- 
teristic time, leading to a revival of the quasihole. 

This dynamics is reminiscent of the collapse and revival 
of a coherent light field which resonantly interacts with a 
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two-level atom. This effect has been studied theoretically 
in the framework of the Jaynes-Cummings model since 
the early 1980s [13, Uli, and has experimentally been ob- 
served in systems of Rydberg atoms [Tol - ITsj , or trapped 
ions [l^ . With the realization of a Bose-Einstein conden- 
sate (BEC) in 1995, also interacting many-body systems 
have become candidates for studying such collapse-and- 
revival effects: In Ref. [20j it has been argued that quan- 
tum fluctuations cause a phase diffusion which leads to 
a collapse of the macroscopic wave function. As a con- 
sequence of the discrete nature of the spectrum, periodic 
revivals of the macroscopic wave function have been pre- 
dicted in Refs. [2l|, It has been proposed to pro- 
duce macroscopic entangled states by time-evolving a 
condensed state j23l . [23 |. An interesting scenario has 
been discussed in Refs. [1^ [2^, studying collapse and 
revival of the relative phase between two spatially sepa- 
rate BECs. Measuring phase correlations between many 
BECs which are distributed on an optical lattice has al- 
lowed for observing the collapse and revival of matter 
waves [13]. Recently, the observation of quantum state 
revivals has been proven to provide relevant information 
about the nature of multi-body interactions in a Bose 
condensed atomic cloud [28j . 

Also the collapse and revival which we discuss in this 
paper allows to extract useful information: The effect 
itself not only clearly distinguishes the Laughlin regime 
from denser ones, but also measuring the revival times 
and positions of the quasiholes allows to determine the 
kinetic and interaction contribution to the energy of the 
system. 

Our paper is organized in the following way: After 
introducing the system in Sec. [TTl we study the coher- 
ent quasihole dynamics in the Laughlin state in Sec. IIIII 
This is in contrast to the collapse-and-revival dynamics 
in denser regimes described in Sec. lIVl In Sec.|Vl we draw 
our conclusions. 



II. THE SYSTEM 

We consider a two-dimensional system of bosonic 
atoms with mass M, described by the effective Hamil- 
tonian H = "Y^^^i^i + V, where the single-particle con- 
tribution reads 

H^ = ' + Y^c'ff G^' + v^)- (1) 

Here, Ai denotes the artificial gauge potential acting on 
the ith particle. It shall describe a gauge field of strength 
B perpendicular to the system. We choose the symmetric 
gauge, Ai = ^[yi,—Xi,Q). Different proposals for syn - 
thesizing this gauge potential are reviewed in Refs. (7l.[l0|. 
The trapping potential is affected by the generation of 
the gauge potentials, but it is possible to make the effec- 
tive trap axial-symmetric with trapping frequency Woff. 

It is useful to introduce a quantity ijJ±_ = \ h^l^ + jfjr, 



which in the case of a rotation-induced gauge field equals 
the applied trapping frequency. From now on it will be 
used to fix units of energy, hw^, and units of length 

The first term in Eq. ([T]) is seen to give rise to Lan- 
dau levels. Using the dimensionless parameter rj = 
B /{2Muj±) < 1, we can express the Landau level gap 
as All = 2?^. The degeneracy of states in each level is 
split by the second term in Hi. In the LLL, the eigenen- 
ergies are given hy Ei = {I — ?])£ + eo, corresponding to 
the Fock-Darwin (FD) states 0f°(z) cx z'^ exp(-|z|2/2), 
with z ~ x + iy, and i the single-particle angular momen- 
tum. The term eo describes an ^-independent zero-point 
energy. The interaction V is assumed to be repulsive s- 
wave scattering, described by 

V=§E^^^H^.-.), (2) 

i<j 

where g parametrizes the interaction strength such that 
gN/'iiT equals the mean-field interaction energy per par- 
ticle. To avoid populating higher Landau levels, we need 
to fulfill gN/Sir <C 77 < 1. Within this constraint, we 
can freely tune either g via Feshbach resonances, or rj by 
modifying the gauge field strength. This allows to drive 
the many-body ground state from a condensate 

N 

vl/o(zi,...,z^)cxn'/'r(^0 (3) 

for <C 1; to the Laughlin state 

^-lIzi, . . . , zat) cx J|(z, - z,f exp(- ^ |z,| V2) (4) 

i<j i 

for 7/ ^ 1 [13. 

After coohng the system into the ground state of H, 
we use laser beams to pierce quasiholes into the state. 
In Refs. [1, [^, it has been confirmed by exact diago- 
nalization that a laser potential focused at position ^, 
Vi oc S^'^^Zi—S,), is able to produce quasihole excita- 
tions in the Laughlin state. The laser intensity / needs to 
be strong enough to close the gap protecting the Laughlin 
state. Defining a quasihole operator Oqh(0 =-^Y\i{^i ~ 
^) which pushes away all particles from position ^ and 
where Af rc-normalizcs the state, the quasihole state cor- 
responding to a state ^'(zi, . . . , z^v) can generally be de- 
fined as 'Sqh{£.;zi,...,ZN) = 6qh(0*(2i- ■ • , -zw)- This 
state has up to N units of angular momentum more than 
the original state. 

If the system is large enough, our procedure can be 
repeated in order to create several distinct quasiholes. 
In the following we investigate the scenario where, after 
generating the quasiholes, the laser beams are abruptly 
switched off. In general, a dynamical evolution is ex- 
pected, since is not an eigenstate of H. This is di- 
rectly clear for ^ 7^ 0, where the operator Oqh(0 breaks 
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the cylindrical symmetry of the Hamiltonian %. We re- 
strict our study to quasiholes characterized by a vanish- 
ing wave function at ^, thus we do not consider the pos- 
sibility of quasiholes of a different form, as for instance 
the half-flux excitations in the Moore- Read state 12911. 



III. COHERENT QUASIHOLE DYNAMICS IN 
THE LAUGHLIN STATE 



First, we consider the Laughlin state with one quasi- 
hole at position ^: 

N 
fc=0 

(5) 

where the fk are totally symmetric polynomials of fcth 
order in the coordinates zi,...,ZAr, with the property 
that each of the coordinates appears at most to lin- 
ear order. Since V^l = 0, we also have V/fc4'L = 0. 
Furthermore, all /jc^'l are homogeneous polynomials in 
the Zi times the overall Gaussian, and thus are eigen- 
states of the single-particle part ^ - Hi with eigenvalue 
Ek = (1 - 77)Lfe + TVeo = (1 - ri)[k + N{N - I)] + Ne^. 
Defining 7 = N{N - 1)(1 - 77) + iVeo and e = (1 - ij), we 
can write the time evolution of the quasihole state: 

N 
fc=0 

(6) 

The exponential in the sum can be absorbed by making 
the ^'s time-dependent: ^(t) = ^c^"s^^*. With this we 
obtain: 

N 



k=0 

= ci(^+^')*vi/L,qh(^>)). (7) 



The overall phase evolution, e'('^+^^)*, is just the dynam- 
ics of a state with a quasihole in the center, ^ = 0, which 
is an eigenstate of As Eq. ([7]) shows, for a symmetry- 
breaking quasihole off the center, the conservation of an- 
gular momentum is ensured by a constant coherent ro- 
tation of the quasihole around origin. Any change in |^| 
is forbidden by conservation of angular momentum. The 
angular velocity of the quasihole is given by — 77). 

Remarkably, no decoherence between different terms in 
the sum of Eq. ([S]) occurs during the evolution. This phe- 
nomenon is a consequence of the state's zero interaction 
energy in a contact potential. Note that this behavior is 
in contrast to conventional fractional quantum Hall sys- 
tems with long-range interactions where a tunneling of 
the quasihole to the edge is expected jsot . 

The above calculation can easily be repeated for more 
than one quasihole. It can generally be shown that 
the time dependence of the corresponding wave func- 
tion can be absorbed into the positions of the quasiholes 



£.1 and an overall phase factor. Therefore 

we just have to note that the normalization factor of 
the wave functions, which carries the information about 
the anyonic statistics of the quasiholes, depends only on 
absolute values and — thus the substitution 
^iS"'^'^* poses no problem there. 



IV. COLLAPSE AND REVIVAL OF THE 
QUASIHOLE 

In contrast to the coherent dynamics of quasiholes 
in the Laughlin state, in this section we will encounter 
coUapse-and-revival processes of quasiholes pierced into 
phases of less angular momentum, L < N{N — 1). 



A. Quasihole in the L — condensate 

We start analyzing the dynamical behavior of a quasi- 
hole in a L = condensate, described by the wave func- 
tion vE'o,qh(C) = Oqh(0*o- As before in the Laughlin 
case, we can decompose this expression into a sum over 
homogeneous polynomials, *o,qh(C) « Z^^^o /^^^"''^o, 
where every term is an eigenstate of the single-particle 
part of T-L, with corresponding eigenvalue fee -|- 7'. Now, 
7' = iVeo is the zero-point energy. We should bear in 
mind that both e and eo depend on 77, so their numerical 
values might be different to the Laughlin case. 

Again, we can absorb the single-particle contribution 
being linear in k into the time evolution of the quasihole, 
so for the non-interacting system, g = 0, we would have 



*o,qh(^,i) = exp 



*o,qh(^~(0) , (8) 



in full analogy with Eq. ([7|). Interactions, however, 
change the situation: The terms (fk = fk'^o are in gen- 
eral not eigenstates of the interaction V. To describe the 
time evolution of this system, we thus have to decom- 
pose the (fk^s into an eigenbasis of V. Since V conserves 
angular momentum, we can restrict ourselves, for every 
(pk, to the subspacc with L = fc, for which we obtain 
the eigenbasis via exact diagonalization. We denote this 
basis by Xa'' and write: 



</5fc 



(9) 



The coefficients cit^ can easily be obtained: The exact 

diagonalization yields the xi'^'' in the Fock basis of occu- 
pation number states \ni^o,nf—i, . . .). In this basis, the 
state (fik is represented by the vector \N — k,k,0, . . .), 
from which it differs only by a normalization factor 
Nk = [{N - k + 1)! {k + We thus have 



Nk 



Ak) 

A" 



N -k,k,0,. 



The Xa^\ being ho- 



mogeneous polynomials of fcth degree, are eigenstates of 
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the kinetic term with the eigenvalue /ce + 7'. We now 
may write 

N 

A:=0 a 

(10) 

Here, ei'^'' is the eigenvalue of V corresponding to the 
eigenvector Xa'^ ■ The presence of this term causes, in 
general, a dephasing of the different contributions to 
'l'o.qh(C)- Thus, while the single-particle contribution 
just rotates the quasihole at fixed radial position |^|, 
the interaction makes the quasihole fade out, as shown 
in the density plots of Fig. [TJ Also, slight deforma- 
tions of the cloud as a whole become apparent during 
the time evolution, in clear contrast to the Laughlin 
case. This interaction-driven dynamics also happens for 
a symmetry-conserving quasihole in the center. 

To quantify the dephasing we consider the eigen- 
values r of the one-body density matrix pij = 
(^o,qh(^7 01 c^Cj IvPo^qhlC: ^)): that is, the occupation of 
the different eigenmodes. Here, Ci (c|) is the annihila- 
tion (creation) operator of a particle in FD state 0f^. 
As shown in Fig. [51 initially all particles occupy the same 
single-particle state, {z — ^) exp(— |2;p/2), and all eigen- 
values of p are zero except for one which is equal to the 
number of particles N. If the hole is placed in the cen- 
ter, the different instantaneous eigenmodes of p are the 
Fock-Darwin functions, as a matter of fact that the initial 
state has definite angular momentum, and the Hamilto- 
nian conserves angular momentum. For all particle num- 
bers we have studied, the most occupied mode is, at any 
time, 01^, while the second most occupied is the Gaus- 
sian, i.e. (j>Q^ ■ As can be inferred from Fig. [21 and as we 
have verified explicitly by calculating the pair-correlation 
function [3l| , the time evolution brings the initially fully 
uncorrelated system into a correlated state, where de- 
tection of one particle at position z influences the out- 
come of another measurement at position z' . The effect 
is most pronounced in small systems with a quasihole in 



the center. In this case, we find for N = 12 a.t t ^ T/2 
a correlated state with almost equally populated modes, 
n^^ = 4.7 and <d = 4.3. 

As seen in Figs. [Hand [H the collapse of the quasihole 
is followed by a perfect revival. We thus have oscillations 
between fully condensed quasihole states and correlated 
states. In the absence of perturbations, these oscilla- 
tions of period T will continue forever. Note that, in 
general, the position of the quasihole at the nth revival, 
changes from period to period due to the kinetic 
contribution. In the following, we will derive expressions 
for T and which will allow to deduce information 
about different system's parameter by measuring these 
quantities. 

Therefore, we have to analyze the spectra which 
define the periods T^''^ ~ 2iTh/e^a^ of the phase oscilla- 

tions of each contribution Xa ■ Now assume that there 
are some energy units u and u', which allow to write 

+ (11) 

with Ha G N. The part linear in k does not depend 
on a, and thus can be absorbed into the kinetic energy 
contribution. It exclusively affects the revival position by 
defining the time dependent parameter ^(t), now rotating 
with an angular velocity (e -I- u')lh. The phase decohcr- 
ence between different a states is controlled by u, and 
the corresponding periods read T^^^ = 2iTh{na^u). We 
see that T = 2Trh/u is a multiple of all Ti''\ so at time 
T, all contributions will have the original phase relations. 
To determine u and w', we numerically analyze the 

(k) 

spectra . We find that the gap above the L — N 
subspace provides us, for any N, with an energy unit 
u « 0.040gA^. All states in the spectrum arc found to 
be given as integer multiples of u plus the ground-state 
energy. Also in subspaces L < N, the same unit u can 
be used to quantize most of the energies. Strikingly, the 
eigenstates to energies which cannot be constructed ac- 
cording to Eq. (jlip have zero overlap with tpk- Thus, 
they do not contribute in Eq. (jTU]). 
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The second energy unit of Eq. turns out to have 
exactly the same value, u' = u. As exact solutions are 
known for the ground state energies in subspaces with 
L < iV 0, |32| . [ssj , we can write down an analytic expres- 
sion u' = gN/Sir. We thus obtain for the revival period 
[in units of wT^l: 



gN 



(12) 



From this formula, we directly see that choosing gN = 
constant makes the oscillation periods independent from 
the size of the system. This choice is convenient as it 
also guarantees a finite interaction energy per particle in 
the thermodynamic limit. Fixing g rather than gN, the 
periods would decrease in larger systems. By Eq. ([T2|) . 
a measurement of the revival period T directly yields 
information about gN. Measuring then the polar angle 
(Z)*^"^ of the revival position ^'^"^ = £,{nT) will allow to 
extract 77. We find 



.) ^ 4n(l ~ ri) 
gN 



(13) 



The effect of the system size can be seen by com- 
paring TV = 8 and iV = 20 at i = T/2 in Figs [1] 
and [D The larger the system, the more it tends to 
maintain its initial properties. In Refs. [1, [33 - l33 |. it 
has been shown that the ground state wave functions 
ioi L < N are closely related to the functions fk, which 
are the wave functions where the total angular momen- 
tum L ~ k < N is most equally distributed amongst 
N particles, fk = J2i<ii<...<i^<N ' " ^^k- From these 
functions, we obtain the polynomial part of the ground 
state wave function of i = fc by replacing the coordinates 
Zi by the relative coordinates Zi = Zi — Z, where Z is the 
center-of-mass coordinate. Note that Z is not just a num- 
ber, but an operator with (Z) = 0. As center-of-mass 
fluctuations decrease with increasing particle number of 
the system, for large-sized systems, Z becomes pinned 
to the center, and the states fk'^o become eigenstates of 
the Hamiltonian with eigenvalues oc k. Just like in 
the Laughlin case, we will then no more observe the col- 
lapse and revival of the hole. The rotational movement 
around the origin will survive the thermodynamic limit if 
the hole is initially placed outside the center. Therefore, 
the dynamics of a single hole does not qualitatively dis- 
tinguish the condensed phase from the Laughlin phase in 
the thermodynamic limit. There is, however, a quantita- 
tive difference, as the period of rotation will be shorter 
than the period of Laughlin quasiholes due to the en- 
ergy u' from Eq. pT|) which has to be absorbed in the 
definition of £^{t). 



B. Piercing two quasiholes in the condensate 

The situation becomes quite different if we pierce a sec- 
ond hole into the condensate. For simplicity, we choose to 
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FIG. 2. Eigenvalues r of tlie one-body density matrix p as a 
function of time. Initially, one hole was pierced at ^ = 0.2 
into the L = condensate, (a) iV = 8. (h) N ^ 20. 



introduce one of them in the center. Our wave function 
is then a linear combination of states \0, N — k, k,0, ■ ■ ■) 
with < fc < A^, and we are in the angular momen- 
tum regime with N < L < 2N. Here, even the GS 
energies El do not behave linearly with L [35|. Also, 
many more states are involved when expressing the states 
\0, N — k,k,0, ■ ■ ■) in terms of eigenstates of the interac- 
tion, which now live in a significantly increased Hilbert 
space. The consequence of this is that, after a quick de- 
phasing, the holes will never exhibit a full revival, see 
Fig. 131 Contrarily, as the comparison of the data for 
N — 8 and A^ = 12 in Fig. |3] suggests, the peaks at the 
same period as given by Eq. are expected to fade 
away quickly in the thermodynamic limit. 

An L = condensate with two quasiholes is similar 
to an L = A" condensate with one quasihole, and in 
the thermodynamic limit the L = N condensate with 
one vortex in the center is the ground state of the sys- 
tem at this angular momentum. Fig. [3] therefore sug- 
gests that in regimes L > N, also a single quasihole 
will dephase. This seems to be reasonable, as sponta- 
neous symmetry breaking has been predicted for states 
with L ^ N [3^ [s^l, and ground states in several 
L subspaces become quasi-degenerate. That means that 
energy differences within each quasi-degenerate manifold 
are very small, while large energy jumps occur between 
different manifolds. Therefore, a dephasing of different 
contributions and thus a collapse of the quasihole must 
be expected. Furthermore, the smallness of energy con- 
tributions in the quasi-degenerate manifolds make revival 
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FIG. 3. Occupation numbers r/N of eigenmodes as a function 
of time in a condensate of = 8 (dashed lines) and A*' = 12 
(solid lines) particles. Initially one hole has been placed in 
the center and second hole at |^| = 1. No full revivals are 
observed, and the peaks at t = nT are washed out in larger 
system. 



times unobservably long. We do not further investigate 
this situation oi L ^ N, as the symmetry-breaking leads 
to the formation of vortex lattices Q , as observed in ex- 
periments . Characteristics of the lattice should allow 
for a clear identification of these phases. 
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FIG. 4. Density profile {x = 0) for a system of A'^ = 6 particles 
initially prepared in the state Oqh(0)^L,qp at different times 
t (in units of cjJ^) with g = \. In (a) the evolution from a 
quasihole in the center to a density peak is shown. In (b) we 
plot the density for all times t < 200, at which it takes an 
instantaneous extremum in the center. 



C. Dynamics of a quasihole in the Laughlin 
quasiparticle state 

Upon increasing the gauge field strength, the vortex 
lattice has been predicted to melt when reaching filling 
factors ^ 6 Q. Then, a variety of strongly correlated 
quantum liquid phases are candidates for the ground 
state. Finally, for = 1/2 or L = N{N-1), the Laughlin 
state becomes the ground state of the system. It is cer- 
tainly in this regime where observable properties to dis- 
tinguish between the phases become most relevant. Let 
us therefore study the dynamics of a quasihole pierced in 
the last incompressible phase which has been predicted 
to occur before reaching the Laughlin state. It is charac- 
terized by L = N(N — 2) and by a wave function which 
differs from the Laughlin wave function only locally at 
the origin, having the form of a Laughlin quasiparticle ex- 
citation, ^-L.qp = exp(- J2i kiP/2)9^i ...dzj^ I\t<jizi - 



For simplicity, we will pierce the quasihole in the ori- 
gin, which makes the resulting state an eigenstate of the 
single-particle part of and all dynamics will exclu- 
sively be driven by V. To obtain the state Oqh,o^L.qp in 



the Fock basis, we numerically diagonalize the Hamilto- 
nian H' = V + Vi in the subspace L = N{N - 1). We 
then decompose this state in the corresponding eigenba- 
sis of V, also obtained by exact diagonalization. Several 
eigenstates of V will contribute, but the largest contri- 
bution comes from the Laughlin state with an overlap 
of 0.709 (0.717) for iV = 7 (iV = 6). Expressed in the 
eigenbasis of V, we can easily perform the time evolution 

of 6qh(0)*L.qp. 

The dynamics is clearly visible in the density shown 
in Fig. |4^: In the course of time, the hole fades out, as 
the center of the cloud gains a finite density. At some 
point, even a density maximum is developed at the ori- 
gin, surrounded by a circular density valley. As the val- 
ley spreads out, the maximum becomes clearly peaked. 
The process is then reversed, and a hole at the center 
re-appears. Such oscillations between a density maxi- 
mum and a density minimum in the center can be ob- 
served repeatedly. The scenario, however, differs from 
the collapse-and-revival process in the condensate: First, 
re-appearing holes are not equivalent to the original hole, 
as their density at the center remains finite, and their core 
size has decreased, see Fig. HJd. Second, the "revival" pe- 
riods are not sharp. In Fig. [Jb, we have chosen precisely 
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FIG. 5. Overlap c; between Oqh(0)4'L,qp and eigenstates of V 
with energy Ei for different A^. 



those times at which the process is reversed. For iV = 6, 
the reversal after the first re-appearance of the quasihole 
is found at i = 55, while the second revival takes place 
at t = 114. 

To understand this behavior, we have to analyze the 
spectrum of V at L = N{N - 2) + TV = A^(A^ - 1). It can 
be divided into a quasi-continuous excitation band and 
the Laughlin state. A gap A separates these two con- 
tributions. For > 6, the gap approaches the constant 
value of A w O.lo, if we choose g (rather than gN) to 
be constant [^,[331 ■ Compared to this value, the energy 
differences between states within the excited band are 
typically very small. This property of the spectrum can 
be seen in Fig.[5l where we have plotted the overlap of the 
eigenstates with the initial state, a = (Oqh(0)'I'L,qp|xi)j 
versus Ei/g. Here, \xi) denotes the eigenstates of V in 
the L = N[N — 1) subspace, and Ei is the corresponding 
eigenenergy. 

Due to this structure, the relative dephasing of differ- 
ent contributions to Oqh(0)^'L,qp from the excited band 
is slow compared to the dephasing of these contributions 
with respect to the contribution from the Laughlin state. 
Following this reasoning, 27r/A k, 63 for g = 1 sets 
the rough time scale for a "quasi" -revival, at which the 
Laughlin state is again "in phase" with the low-energy 
contributions from the excited band. This slightly dif- 
fers from the number we find by analyzing the density 
(cf. Fig. Hb), t = 55 for = 6. But we note that the 
most important contribution from the excited band is a 
state with energy E = 0.112,9 (cf. Fig[S]). Thus, it is 
in phase with the Laughlin after a time t ~ 56. It has 
an overlap with Oqh(0)^L.qp of 0.533. A superposition 
of this state and the Laughlin state is able to reproduce 
the quasihole state with a fidelity of 78%. Other impor- 
tant states, with overlaps 0.279 and 0.110, are found at 
E = O.lOSg and E = 0.102.g. At i = 56 they are stih 
nearly in phase with the E — 0.112^ state. But also 
states with E = 0.1645 and E — 0.182^ contribute sig- 
nificantly with overlaps 0.180 and 0.133. These states 



will be clearly out of phase, making the revival imper- 
fect. Subsequent revivals will more and more suffer from 
the slow dephasing within the manifold of excited states. 
This explains the small irregularity in the revival peri- 
ods and the loss of the quasihole character in the density 
profiles (see Fig. |3Jd). 

Posing the question whether the described dynamics 
will survive in the thermodynamic, we first note that 
due to the similarity between the Laughlin state and 
Oqh(0)^L.qp, we can always expect contributions from 
both the Laughlin state and the excited band. This as- 
sessment seems to agree with Fig. [51 where the overlap 
with the Laughlin state is found to be almost constant 
while varying particle number, 5 < < 7. As the 
Laughlin gap is known to be constant for large A^, also 
the revival time has to be. The imperfection of the re- 
vival, characterized by a finite density at the center, con- 
tinuously improves, as we increase the system size from 

= 4 to A^ = 7. This seems reasonable, since energy 
differences in the excited band decrease with larger A^, 
slowing down the dephasing in the excited band. In Fig. [5] 
this reflects in the decreased spreading of relevant states 
when increasing particle number. However, whether this 
might lead to a perfect collapse-and-revival process in the 
thermodynamic limit, is not clear from our calculations. 

V. CONCLUSIONS 

We have shown that observing the dynamics of a quasi- 
hole might serve to classify different ground states in 
the LLL regime. Especially, the absence of decoher- 
ence is a characteristic feature of the Laughlin regime, 
L > N{N — 1), due to its zero interaction energy. It 
is in clear contrast to the collapse of the quasihole ob- 
served in finite systems with L < N{N — 1). For sim- 
plicity, we have considered an idealized system with a 
cylindrically-symmetric Hamiltonian. We note however 
that even deformed Lau ghli n states, as expected in laser- 
induced gauge fields (ol. Il2|. should be characterized by 
a decoherence-free quasihole dynamics due to their van- 
ishing interaction energy. In the condensed regime, the 
collapse of the quasihole is followed by a perfect revival, 
and the system oscillates between a condensed and a cor- 
related state. System parameters like gN and rj, speci- 
fying the interaction and the single-particle energy, can 
be obtained by measuring period and positions of this re- 
vival. A collapse-and-revival is also found for a quasihole 
in the Laughlin-quasiparticle state, being an incompress- 
ible phase in the direct vicinity of the Laughlin phase. 
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